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ABSTRACT 

Type IIP (Plateau) Supernovae are the most commonly observed variety of core collapse events. 
They have been detected in a wide range of wavelengths from radio, through optical to X-rays. 
The standard picture of a type IIP supernova has the blastwave interacting with the progenitor's 
circumstellar matter to produce a hot region bounded by a forward and a reverse shock. This region 
is thought to be responsible for most of the X-ray and radio emission from these objects. Yet the origin 
of X-rays from these supernovae is not well understood quantitatively. The relative contributions of 
particle acceleration and magnetic field amplification in generating the X-ray and radio emission need 
to be determined. In this work we analyze archival Chandra observations of SN 2004dj, the nearest 
supernova since SN 1987A, along with published radio and optical information. We determine the 
pre-explosion mass loss rate, blastwave velocity, electron acceleration and magnetic field amplification 
efficiencies. We find that a greater fraction of the thermal energy goes into accelerating electrons than 
into amplifying magnetic fields. We conclude that the X-ray emission arises out of a combination of 
inverse Compton scattering by non-thermal electrons accelerated in the forward shock and thermal 
emission from supernova ejecta heated by the reverse shock. 

Subject headings: Stars: Mass Loss — Supernovae: Individual: SN 2004dj — shock waves — circum- 
stellar matter — radio continuum: general — X-rays: general 



1. INTRODUCTION 

Core collapse supernovae with hydrogen lines near 
maximum light and a pronounced plateau in their vis- 
ible band light curve that remain within ~ 1 mag of 
maximum brightness for an extended period are called 
type IIP supernovae. The plateau duration is often 60- 
100 rest-frame days and is followed by an exponential 
tail powered by radioactive decay at late times. Type 
IIP supernovae constitute about 67% of all core collapse 
supernovae in a vo lume limited sample (d < 30 Mpc) 
(|Smartt et al.l I2009T ) . Their characteristic optical light 
curves are attributed to the hydrogen envelope of the 
progenitor remaining largely intact before the core col- 
lapse. Several lines of evidence suggest that these stars 
were red supcrgiants when they exploded. Type II super- 
novae (as well as the type Ib/c supernovae) are associ- 
ate d with recent star forming regions of the spiral galax- 
ies dFilippc nko 1997]) which suggest that their progenitor 
stars are relatively short lived and therefore explode from 
massive stars (M > 8Mq) which have faster burning of 
their nuclear fuels at generally higher central temper- 
atures compared to less massive stars. Their plateau 
brightness and duration combined with their expansion 
velocities suggest a pre-supernova radius of typically red 
supcrgiant dimensions (e.g. 10 2 ~ 3 Rq). While there 
are more massive red supcrgiants in the Local Group of 
galaxies, no high-mass red supcrgiant progenitors above 
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17M have been found for IIP supernovae in direct detec- 
tion efforts in volume-limited sea rch for supernova pro- 
genitor stars (jSmartt et al.l I2009T) . This red supergiant 
problem has led to a viewpoint that these more massive 
stars could have core masses high enough to form black 
holes and supernovae from them could be too faint to 
have been detected. At the same time a minimum stellar 
mass for a SN IIP to arise from is about 8.5 ±1.0 Mq, 
consistent with the upper limit to white dwarf progenitor 
masses. A census of type IIP progenitors and the par- 
tial loss of their hydrogen envelopes could potentially be 
important for neutron star vs black hole formation and 
their number distribution in the galaxy. 

In a type IIP supernova, the expanding ejecta interacts 
with the slow wind of the red supergiant even though the 
star retained most of its hydrogen envelope intact. This 
interaction generates a hot region, bounded by the for- 
ward and reverse shocks, which may emit thermal X-rays. 
These shocks may also accelerate charged particles and 
synchrotron radiation from relativistic electrons can lead 
to radio emission that is strong enough to be detected at 
cxtragalactic distance scales. The measurable radio and 
X-ray properties give complementary information about 
the regions shocked by the forward and reverse shocks 
and can constrain the physical properties of the interac- 
tion region and the progenitor star. Because type IIP 
supernovae have a plateau phase of high optical lumi- 
nosity, energetic electrons near the forward shock find a 
dense photon environment of seed photons which may 
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be Compton boosted to X-ray energies. The resultant 
Compton cooling can affect the population of the elec- 
trons which also emit at radio wavelengths and show sup- 
pressed radio fluxes above a characteristic cooling break 
until the end of the plateau phase where the optical pho- 
ton d ensity undergoes a rapid decline ([Chevalie r et al.l 
12001 . 

Using SN 2004dj as a prototype of its class we would 
like to address the questions: 

1. What is the origin of X-rays detected by the Chan- 
dra from young type IIP supernovae? 

2. How fast was the progenitor losing mass via winds 
in the final phase before explosion? 

3. How efficiently does the supernova accelerate cos- 
mic ray electrons in its forward shock? 

4. What is the extent of turbulent magnetic field am- 
plification in the post-shock material? 

5. Are synchrotron and inverse Compton losses im- 
portant in explaining the radio lightcurves? 



SN 2004dj, discovered by K. Itagaki (see lNakano et al.l 
12001 on 31.76 July 20 04 (UT) in the spiral gala xy NGC 
2403, was the closest (iKarachentsev et aT1l2004f) normal 
type IIP SN (jPatat et al.l I2004D observed to da te. It 
was discovered 50 ± 21 days ( Zhang et al.l 12006^ 1 after 
its progenitor's core collapse and had a peak m agnitude 
of 11.2 mag in the V-band (|Zhang et all 120061) . Anal- 
ysis of its pre-explosion images indicated that the ini- 
tial main sequence mass of the progenitor of SN 2004dj 
was about 15 Mp, dMafz-Apellaniz et aT1l2004[ ): see how- 
ever lWang eTaTTpOOl " where the progenitor had a main- 
sequence mass of 12M Q . SN 2004dj occurred at a posi- 
tion coincident with Sandage star number 96 (hereafter 
S96, and resolved by HST to be a young star cluster) in 
NGC 2403. To date, only seven supernovae of type IIP 
are known X-ray emitters. Only four of them are known 
radio emitters: Supernovae 1999em, 2002hh, 2004dj and 
2004ct. X-rays from circumstellar interaction has been 
studied i n detail for type II P supernovae i ncluding SN 
1999e m (jPoolev et al J 12002ft . and 2004et (jMisra et al.l 
|2007|) . Because of its close range, SN 2004dj was eas- 
ily detectable in the X-ray and radio. The Chandra X- 
ray Observatory had targeted it multiple times early af- 
ter its explosion . The e arliest detection was reported by 
iPoolev fc Lewinl (|200l . 

The X-ray and radio luminosities arc both sensitive 
to the mass loss rate and the initial mass of the pro- 
genitor star. In this paper we demonstrate that for SN 
2004dj the determination of the inverse Compton pow- 
erlaw spectrum from the forward shock, the characteris- 
tics of the X-ray line emitting thermal plasma from the 
reverse shock and the radio light curve peak, uniquely 
determines some properties of the explosion and the 
progenitor star. These include the pre-supernova mass 
loss rate (M), the fraction of post-shock energy density 
which goes into relativistic electrons (e e ) magnetic fields 
(es). The time dependent nature of the thermal and 
non-thermal X-ray fluxes also bear signatures of the cir- 
cumstellar density profile as the forward shock encoun- 
ters more and more circumstellar matter and the ejecta 
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Fig. 1. — First X-ray spectrum of SN 2004dj on 2004 August 
9. Bars are counts from Chandra, binned for display. Dotted line 
is the powerlaw model for the inverse Compton flux. Dashed line 
is the APEC model for the thermal plasma. Solid line is the full 
model. Note that the full model is dominated by the non-thermal 
flux even at ~ 1 keV. 



TABLE 1 

Observation summary of SN 2004dj with Chandra 



Date 


Exposure 


Count Rate 


Flux (0.5-8 keV) 


(2004) 


(ks) 


(10~ 3 sec- 1 ) 


(ergs cm" 2 s — 1 ) 


Aug 09 


40.9 


12.80 ±0.56 


8.81 X 10~ 14 


Aug 23 


46.5 


10.03 ±0.47 


6.98 x 10~ 14 


Oct 03 


44.5 


5.60 ±0.36 


3.30 x 10- 14 


Dec 22 


49.8 


3.05 ±0.25 


2.02 x 10~ 14 



profile as the reverse shock ploughs into the expanding 
ejecta. We also see the correlation of the inverse Comp- 
ton flux with the optical light curve. 

2. OBSERVATIONS OF SN 2004DJ 

As the nearest (~ 3 Mpc) supernova since SN 1987A, 
SN 2004dj enjoyed extensive multi-wavelength coverage. 
The details of the observations in Radio, Optical and 
X-rays, used in this work, are given below. 

2.1. Chandra X-Ray Observations 

SN 2004dj was observed by a Target of Opportunity 
program (PI: Walter Lewin, Cycle: 5, ObsIDs: 4627- 
4630) using the Chandra X-ray Observatory on four oc- 
casions: 2004 August 9 (See FigHJ, August 23, Octo- 
ber 3 and December 22 (See Fig©. The ACIS-S chips 
were used on all occasions, without any grating, for 50 
ks each. The supernova was detecte d in the first of these 
observations (|Poolev fc L ewin 2004]) . We analyze all four 
observations in this work for the first time. See Table [1] 
for details of observations. 

Before spectral analysis, the data from separate epochs 
were processed separately but in a similar fashion, fol- 
lowing the prescribed threads from the Chandra Science 
Center using CIAO 4.4 with CALDB 4.4.8. The level 
2 events were filtered in energy to reject all counts be- 
low 0.3 keV and above 10 keV. The resulting events were 
mapped and the supernova easily identified. The source 
region was marked off and a light curve was generated 
from the rest of the counts. Flare times were identified in 
this light curve, as time-ranges where the count rate rose 
above 3 times the rms. This resulted in a good time inter- 
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val table which was used to further filter out the events. 
The spectrum, response and background files were then 
generated for these good events. In order to obtain the 
highest available spectral resolution, the data were not 
binned; unbinned data were analyzed in the next steps. 

2.2. Radio Observations 

SN 2004dj has been observed extensively in radio bands 
from Aug 2004 until May 2007 with various telescopes 
including the Very Large Array (VLA), Giant Metre- 
wave Radio Telescope (GMRT) and Multi-Element Ra- 
dio Linked Interferometer Network (MERLIN). MERLIN 
observed SN 2004dj extensively in the 5 GHz band cover- 
ing the period from 2004 Au gust 5 to December 2. These 
observations are reported in iBeswick et a l. (2005). The 
GMRT observations started on 2004 Aug 12 and con- 
tinued until 2007 May 22. The observations were made 
in the 1420, 610, 325 and 235 MHz bands. The first 
t wo GMRT observations in 1420 MHz band are reported 
in lChandra fc Ravi (pOOl . The VLA followed SN 2004dj 
extensively starting from Aug 1, 2004 until May 28, 2007. 
The observations were made in all VLA bands starting 
from 1.4 GHz band upto 44 G Hz band. The f i rst 8 4 
GHz observation was r eported bv lStockdale et all (|2004D . 
IChevalier et al.l (|2006h discussed these published radio 
observations of SN 2004dj to interpret the physical prop- 
erties of SN 2004dj and its interaction with the circum- 
stellar matter. They established Free Free Absorption 
to be the radio absorption mechanism based on the pub- 
lished radio observations. A comprehensive paper in- 
cluding all the radio observations and their detailed in- 
terpretation is underway (Chandra et al., 2012, under 
preparation.) 

3. A TALE OF TWO SHOCKS 

Supernova ejecta hits the pre-explosion wind at a ve- 
locity much larger than its characteristic sound speed. 
This creates a strong forward shock moving into the 
circumstellar m atter arising fro m the mass loss from 
the progenitor (jChevalierl 119821 ) . In case of a type 
IIP supernova, the material behind this shock is too 
hot (~ 100 keV) and tenuous to contribute signifi- 
cantly to the X-ray flux in the Chandra bands (0.3 — 10 
keV) through its thermal emission. However, the ra- 
dio emission from type IIP supernovae is said to arise 
from the non-thermal elec trons accelerated in this re- 
gion (jChevalier et al.ll2006D . Some optical photons from 
a supernova may also be inverse Compton scattered 
off energetic electrons in this region into the X-rays 
(jSutaria et al.ll2003t iBiBrnsson fc Franssonl 120041 ) . Scat- 
tering from a non-thermal population of electrons with a 
power law index p generates a spectral index of {p — l)/2 
in the optically thin radio synchrotron. So the inverse 
Compton spectrum may be modeled as a power law with 
a photon index of (p + l)/2. 

In the frame of the ejecta however, the circumstellar 
m atter impinges supersonically on the ejecta. According 
to IMcKed (|1974[ ) this drives an inward (in mass coordi- 
nates) propagating wave which reheats the ejecta, which 
would have otherwise been adiabatically cooled by its 
rapid expansion, to around ~ 1 keV. A supernova typi- 
cally ejects more matter at lower velocity. As the reverse 
shock moves into slower moving ejecta, its heats up more 
gas which can contribute to the X-ray flux. The plasma 
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Fig. 2. — Last X-ray spectrum of SN 2004dj on 2004 December 
22. Note that the relative contribution of the non-thermal inverse 
Compton, with respect to the flux from the thermal plasma has 
gone down, when compared with the first epoch. The full model is 
now dominated by the thermal flux at ~ 1 keV. 

in this region can be modeled as a collisionally-ionized 
diffuse hot gas. 

3.1. X-ray Spectral Fitting 

We imported the extracted spectra into XSPEC 12.7.1 
for further analysis. The data from all epochs were 
jointly fitted with the combination of a simple photon 
power law using powerlaw a nd collisionally-ion ized dif- 
fuse gas using APEC models (|Smith et al.ll200l[) , passed 
through a photo-electric absorbing column using wabs. 
The column density for wabs, plasma temperature for 
APEC and photon index for powerlaw were shared be- 
tween all epochs. The emission measure for APEC and 
normalization for powerlaw were determined separately 
for each epoch from the joint analysis. As, the fits were 
performed on unbinned data, each bin would have too 
few photons for the x 2 statistic to be useful . We t here- 
fore used the cstat statistic based on iCashl (|1979f ). In 
order to evaluate the goodness of fit, we simulated 10000 
spectra based on the best fit model and found out that 
only 25% of these simulations had cstat less than that for 
the data. We therefore conclude that this model provides 
a good fit to the data. The best fit absorption column 
density is n H = (1.7 ± 0.5) x 10 21 cm" 2 . See Tablc[T]for 
model fluxes on each date. 

3.2. Reverse-Shocked Material 

A self-similar solution for the interaction of the fast 
moving ejecta and the nearly sta tic (neglig i bly m oving) 
circumstellar matter was found bv lChevalierl (|1982D . This 
can be used to generate a complete description of t he 
forward and reverse-sh o cked m aterial (|Nadezhinlll985f ). 

Following lChevalie"rl (|1982t ) we consider supernova 
ejecta described by 

Pcj CXi- 3 ^", (1) 

where the density falls due to free expansion and higher 
velocities have less matter, described by the power law 
index rj. This ejecta interacts with the circumstellar mat- 
ter, described by a power law profile p cx r~ s . For a 
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Fig. 3. — Confidence contours for Temperature (in the APEC 
model) vs Photon Index (in the powerlaw model). Note that the 
Chandra data can successfully break the degeneracy between the 
two components and account for their most important parameters, 
namely the temperature and photon index. 

steady wind (s = 2), we have 



.4 



M 



t.2 4:TTr 2 v v 



(2) 



where M is the mass loss rate and tj w is the w ind veloc- 
ity o f the progenitor. Here A = M/(47rv w ) (jChevalierl 
1982). Shock jump conditions for a fast shock dictate 
that the density behind the forward shock is given by 
p cs = 4p w . In the thin shell approximation, this can be 
further related to the density behind the reverse shock 

fa-3)(77-4) 



Po 



(3- S )(4- S )< 



(3) 



According to iMatzner fc McKed (|1999t) 77 is expected to 
be 12 for a red supergiant progenitor. Using this value 
simplifies the above expression to po = 36p cs = 144p w . 

Therefore the reverse-shocked material is much denser 
than the circumstellar matter. So it can cool efficiently 
via radiative processes and contribute to the thermal X- 
ray flux. 

3.3. Blastwave Velocity 

The velocity of a shock determines the energy im- 
parted to particles crossing it, and hence the post-shock 
plasma temperatures. Furthermore, the forward and re- 
verse shock velocities are connected by jump conditions 
once the ejecta and circumstellar matter profile are spec- 
ified. So, the post-shock temperature of ejecta which has 
crossed the reverse shock can be expressed in terms of the 
blastwave velocity. Since the former is an observable, the 
la tter can be der i ved fr om it. 

INvmark et al.l (|2006D give the temperature of the 
reverse-shocked material To, as 



To = 2.27 x 10V (3 S) 



(v - s ) 



2 vi K, 



(4) 



where p is the mean atomic weight per particle, r\ is 
the power law index of the ejecta profile, s is the power 
law index of the circumstellar density profile and V4 is 
the forward shock velocity in units of 10 km s _1 . p is 
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Fig. 4. — Confidence contours for normalization of the thermal 
flux (in the APEC model) the normalization of the non-thermal 
flux (in the powerlaw model). Note the anti-correlation between 
the thermal and non thermal fluxes, as their sum must explain the 
total flux which is fixed by observations. 

0.61 in the solar-like outermost zone of the progenitor 
(|Nvmark et al.ll2006D . s is 2 for a circumstellar matter 
set up by a steady wind, and we again use T] = 12. With 
these assumptions, the forward shock velocity can be ex- 
pressed as 



V s = 10 4 



fcT 



kms 



(5) 



1.19 keV 

Since, the best fit temperature is 0.997 ±0.054 keV (See 
Fig©, the implied velocity is V s = (9.2 ± 0.3) x 10 3 km 
s _1 . 

This is consistent with the fastest-moving ejecta in a 
typical type IIP supernova. Therefore, the observed ther- 
mal component of the X-ray spectrum can be explained 
as emission from reverse- shocked mater i al, pr ovided this 
is the blastwave velocity. iChugai et al.l (|2007[ ) have pro- 
posed diagnostics for circumstellar interaction in Type 
IIP supernovae using the detection of high-velocity ab- 
sorption features in Ha and He I lines during the pho- 
tospheric stage. The highest velocity Ha absorption fea- 
ture observed in the optic al spectra of SN 20 04dj is at 
8200 km s" 1 between 6 4 (|Chugai et al.l |2007[) and 102 
days (|Vink6 et al.l I2006T ) wit h respect to t he ex plosion 
date of June 28 proposed bv IChugai et al.l (|2007[ ). This 
overlaps with the epoch of Chandra observations and 
provides a lower limit to the allowed forward shock ve- 
locity, so the velocity that we obtain is consistent with 
constraints from optical spectroscopy. 

3.4. Circumstellar Density 

The reverse-shocked matter contributes bulk of the 
thermal X-rays. In order to find the circumstellar density 
it is necessary to estimate the amount of reverse-shocked 
material in terms of the forward-shocked material. Using 
Equations [5] and [3] and again assuming 77 = 12, we get 

36M 

Pa = — 5 — • (6) 
irr z v,„ 

Similarly, following INvmark et all (2006) we can express 
the total mass of reverse-shocked material as 



Mr, 



-Mr, 



UIR S 



(7) 
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where R s is the blastwave radius. Multiplying these two 
results, we get, 



144 / M \ 1 
poMo = — — . 

7T \ V w R s 



(8) 



These have to be converted to number densities to get 
the emission measure which determines the thermal X- 
ray flux. For a solar-like composition appropriate for 
the outermost shells of the supernova, p — 1.17amu x 
n e = 1.40amu x n#. Using these, we have the emission 
measure for the reverse-shocked material as 

_ 144 / M\ 1 
n e n H dV -—y—j - 17amu)(1 AQamu)Rs (9) 

We estimate the shock radius from the velocity deter- 
mined in the previous Sub-Section as R s = V s t, with 
the time calculated from th e explosion date of 2004 June 
11 ± 21 (jZhang et al.ll200l or JD 2,453,168. As all the 
other quantities are known, the circumstellar density can 
now be determined from the emission measure. 

The APEC fit to the thermal part of the X-ray spec- 
trum of SN 2004dj gives, 

io- 14 



A-k\D a {1 + z)]- 



J n e n„dv = (7.5 ±2.7) x 10~ 6 , (10) 

for the first epoch (see Fig|4}. NGC 2403 is nearby and 
has negligible redshift. We adopt a dis tance to NGC 240 3 
of 3.06 Mpc based on Cepheids from IS aha et al.l (|2006l ). 
Eliminating the emission measure between the last two 
equations, we have, 



— = (2.0 ±0.7) x 10 13 g s" 1 . 
Therefore the mass loss rate of the progenitor was 



M= (3.2 ±1.1) x 10~ 



M yr" 



(11) 



(12) 



10 km/s, 

This gives us A; however IChevalier fc Franssonl (|2006fl 
define a non-dimensional = A/(h x 10 11 g cm" 1 ). For 
SN 2004dj, we therefore get A* = 3.2 ± 1.1, which can 
now be compared with estimates for circumstellar density 
from the radio emission. 

3.5. Forward- Shocked material 

Electrons are believed to be accelerated to relativistic 
velocities in the forward shock. These are in turn, re- 
sponsible for the r adio emiss i on fro m supernovae. The 
simplest model by Chevalier (1982), is to assume that 
a fraction e e or eb of the thermal energy is used to 
accelerate electrons and amplify magnetic fields respec- 
tively. As a result the radio emission from a supernova 
is dependent upon these fractions and does not mea- 
sure t he circumstellar density direc tly. Instead, accord- 
ing to IChevalier fe Franssonl (|2006l) . 



the radio emission 



measures 



5, =A*e B -ia 8/19 = 1.0 



i 



-8/19 



o.5y V mJ yy 



-4/19 



\Mpc) V5GHz 



-4/19 



(13) 



where F p is the peak flux at peak frequency v at 10 x 
iio days from the explosion and e^-i = es/0.1. The 
cquipartition factor is defined as a = t P Jen- Us i ng the 
4.99 GHz, radio light curve from iBeswick et all ()2005l ) 
for SN 2004dj, we have S* = 5.1. 

The same relativistic electrons which generate the ra- 
dio spectrum via synchrotron emission may also con- 
tribute to the X-ray flux by inverse Compton scattering 
the optical photons from the supernova's peak bolomet- 
ric luminosity of around ~ 10 42 ergs s -1 . This may ex- 
plain the dominant non-thermal component seen in the 
X-ray spectrum at early times (see Fig [lj . For an elec- 
tron index p = 3 this is expected to generate a powerlaw 
with photon index 2, consistent with the o b serva tions 
(see Fig [3]). Following IChevalier fc Franssonl (|2006ft the 
normalization of the inverse Compton flux at 1 keV can 
be written as 



E 



dE 



^ m36„, o „H/19t/ 
X 10 7min<J*a Vi 



/ L bol (t) 
\10 42 ergs s _1 



tio ergs s \ 



(14) 



where 7 m i n is the minimum Lorentz factor of the rela- 
tivistic electrons. 

The normalization of the inverse Compton flux ob- 
tained for the first epoch (ii ~ 6) of Chandra observa- 
tions, is seen to be (3.8 ± 0.5) x 10 37 ergs s _1 (see Fig 
Substituting this in the left hand side and using 
a = 5.1 as found in this work, V4 = 0.92 as im- 
plied by the temperature of the reverse-shocked pla sma, 
L hoX = O.i 

a^23x 7mi ^ n , 



where 7 m i n = 1 if the electron spectrum extends all the 
way down to those at rest. However, if one only con- 
siders relativistic elec t rons w ith 7 m j n ~ 2.5 as used by 
IChevalier fc Franssonl (|2006ft for SN 2002ap, it implies 
a ~ 4.8. We prefer this latter (more realistic) assump- 
tion as the bulk of the accelerated electrons cannot be 
at rest or very low velocities. So we shall use this latter 
value in further calculations. 

Now that we know S+ from radio synchrotron, A+ from 
thermal X-rays and a from inverse Compton, we can use 
the definition of S+ to get 



e B = 0.1 x 6 B -i - 0.082, 
which in turn implies 

e e = a x e B ~ 0.39 . 



(16) 



(17) 



Hence we have determined both e e and cb using a com- 
bination of X-ray, radio and optical data. 

4. TIME EVOLUTION 

The first things that are apparent from the Chandra 
observations of SN 2004dj is that the total X-ray flux 
softens and falls with time. These gross features have 
already been observed in anoth er type IIP supernova, 
SN 2004et bv iRho et all (|2007f l. However, the quality 
of the SN 2004dj spectra allows a quantitative expla- 
nation for the first time. While both the thermal and 
non-thermal components fall eventually, their propor- 
tion changes with time. At early times, the spectra is 
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Fig. 5. — Time evolution of the normalization of the thermal 
APEC component (proportional to the emission measure) in the 
spectrum of SN 2004dj obtained from simultaneous fits to Chan- 
dra spectra. The time interval is measured with respect to the 
explosion date of June 11 from [Zhang ct al. ( 2003) . Line is best 
fit ccT 1 function. Note that the data do not rule out a constant 
mass loss rate model, but suggest one with variable mass loss. 

dominated by the inverse Compton flux and is therefore 
harder (see Fig [TJ . As the bolomctric light curve of the 
supernova decays, the source of seed photons to be scat- 
tered turns off. So at late times, the spectra is dominated 
by emission from the reverse-shocked thermal plasma and 
is therefore softer (see Fig [2]). 

4.1. Thermal X-Rays 

Since the post shock temperature in a self similar blast- 
wave does not change much, as the shock velocity is 
slowly varying, the variation in the thermal X-ray flux 
would come from the variation in the emission measure. 
We have already derived in Equation[5]that the emission 
measure scales as 1/R S . Therefore, if the radius scales 
approximately as R s = V s t, then we have 



n e riHdV oc t 



(18) 



for a constant pre-explosion mass loss rate. We compar- 
ing this to the temporal variation of the normalization 
of the APEC flux from SN 2004dj (see FigEJ). While a 
constant mass loss rate cannot be ruled out, there is a 
hint of variation up to 50% from a constant mass loss. 

4.2. Inverse Compton 

The time variation in the inverse Compton flux 
comes from the expansion if the blastwave and the 
dimming of the sup e rnova 's supply of seed photon. 
IChevalier fc Franssonl ([2006D have shown that the nor- 
malization of the inverse Compton flux scales as 

E dL I c iboiW ( } 

dE t y ' 

Thus the inverse Compton flux is expected to behave as 
oc i -1 in the plateau phase where Lbol is nearly constant 
for a type IIP supernova. Thereafter, as the supernova 
luminosity falls off, the inverse Compton flux is expected 
to fall off from this behavior. This prediction is borne 
out by the observations as the first two epochs (see Fig 
[5]) which lie in the plateau phase can be fitted by a oc t -1 
function, while the late time flux normalization falls off 
significantly below this curve. 




60 80 100 120 140 160 180 200 
t (Days since explosion) 

Fig. 6. — Time evolution of the normalization of the non-thermal 
inverse Compton component in the spectrum of SN 2004dj obtained 
from simultaneous fits to Chandra spectra. Line is best fit oc t _1 
function. The time interval is measured similarly as in Fig[5] Note 
that the first two epochs, for which Lbol i s nearly constant, lie on 
the line. Also note that one the L bo i decreases due to the termi- 
nation of the plateau phase beyond 100 days, the inverse Compton 
flux also falls off. 

5. IMPLICATIONS 

Our results have wide ranging implications for the pro- 
genitors of type IIP supernovae, conditions of the plasma 
in the supernova shock waves, particle acceleration and 
magnetic field amplification schemes. These are dis- 
cussed briefly below. 

5.1. Progenitors of type IIP Supernovae 

According to iMaiz-Apellaniz et al.l (|2004j ) the progen- 
itor of SN 2004dj had a mass of ~ 15M©. In this 
mass range the progen itor dies as a red supergiant and 
IChevalier et al.l (|2006l ) find theoretical mass loss rates 
in the range M = (0.84 - 1.6) x 10~ 6 M yr" 1 for 
such progenitors. However, we find pre-explosion mass 
loss rate for the progenitor of SN 2004dj to be M — 
(0.32 ± 0.11) x 10~ 6 (vw/10 km s" 1 ) M Q yr" 1 . This can 
be reconciled if the wind velocity is > 26 km s _1 , which 
is unlikely as the wind velocity of a typical red supergiant 
is around 10 — 15 km s _1 . However, these r ates have been 
calcul ated for stars of solar metallicity. iSmartt et al.l 
(|2009f ) suggest that the metallicity of the SN 2004dj site 
is log[0/H] = 8.4. Assuming a solar value of 8.7, this im- 
plies a h ost metallicity of Z = 0.5Z W . Using M oc Z 5 
following lSchaller et alj (|1992l ); lHeger et al.l (|2003D would 
bring down the theoretical values closer to the observed 
one. 

Stars in the mass-range 1 1-19 Mq were evolved using 
MESA (jPaxton et al.H2011l) for three different metallici- 
ties z=0.4, 0.50 and 0.75 Zq. The mass loss scheme that 
was used is described as the Dutch Scheme in the MESA 
code; see Sec 6.6 MESA code paper by iPaxton et al.1 
(|2011l ). This scheme switches on a RGB wind depend- 
ing upon the burning stage. The scheme takes care of 
the surface temperature changes in the different stages 
of stellar evolution and adjusts itself accordingly. The 
mass loss in the RGB branch, which is rel evant for the fi- 
nal sta ges of type IIP progenitors, follows Ide Jager et al.1 
(|1988f h The stars were evolved till the central density 
had reached 10 12 g cm~ 3 . The calculated mass loss was 
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Fig. 7. — Zero Age Main Sequence Mass (ZAMS) and Mass Loss 
rate (in the last decades of their lives) for the progenitor of SN 
2004dj (Magenta Cross), MESA runs of 0.4 (solid red), 0.5 (dotted 
blue) and 0.75 (dashed green) Zp, metallicity, and theoreti cal line 
(dot-dashed cyan) from Nicuwcnhuiizcn & do Jagcr (1990) (with 
R = 10 3 i?©) plotted for comparison. Error bars are 3<r confidence 
intervals. Note the disagreement between the observed properties 
of the SN 2004dj progenitor and the theoretical models. 

the average value over the last century of evolution. This 
is the mass which the blastwave sweeps up during first 
few years of its evolution. 

In Fig[7]we co mpare the progenitor mass of SN 2004dj 
as determined bv iMaiz- Apellaniz et a l. (2004) and mass 
loss rate determined in this work, with the expectation 
from stellar evolution and standard mass loss rate pre- 
scriptions. There is an obvious disagreement between 
the observed mass loss rate and theoretic al predictions. 
Such a disagreement was also noted by iChugai et al.l 
(|2007t ) who determined mass loss rates for the SN 1999em 
and SN 2004dj using absorption lines from optical spec- 
troscopy. 

5.2. Ionization Equilibrium 

The analysis of X-ray spectra to yield parameters for 
the underlying explosion models of young supcrnovae and 
SNRs is significantly complicated by the possibility that 
their X-ray emitting plasmas may not be ionization equi- 
librium. The X-ray emission comes from an impulsively 
heated (by the blastwave or the reverse shock) uniform 
and homogeneous gas which is initially cold and neu- 
tral. For young extragalactic supernovae, which are point 
sources contained within the telescope beam-size, the X- 
ray spectra are spatially averaged over the whole object. 
In the case of non-equilibrium plasma the spectra depend 
upon the shock velocity V s , the ionization timescale To 
(the product of the supernova's age t and the post-shock 
electron number density n e ) and the extent of electron 
heating at the shock front (which determine how differ- 
ent the electron and ion temperatures T e ,T s are). Spec- 
tra consisting of lines whose strengths depend upon the 
heavy elements abundances and their relative fractions. 
Both X-ray continuum and lines arise from electrons un- 
dergoing bremsstrahlung or elec tron impact excitation 
of ions. iBorkowski et al.l ()2001[ ) have emphasized that 
for SNRs ionization equilibrium models are often not ad- 
equate, and the constant-temperature, single-ionization 
timescale NEI model is better than an equilibrium model. 
Young supernovae have shorter lifespans till their present 



age but may be expanding through denser circumstcllar 
media and so these issues of adequacy of NEI vs IE mod- 
els may or may not be relevant. 

Because substantial X-ray data from Chandra exists 
in multiple epochs for this target, we investigate the 
relevant parameter space to see whether NEI models 
would constrain emission model parameters sufficiently 
well to require their use (as opposed to IE models). 
The NEI model is an XSPEC provides a non-equilibrium 
ionization collision al plasma model which uses APED 
([Smith et al.ll200ll ) to calculate the resulting spectra. It 
assumes a constant temperature and single ionization 
timescale parameter. In NEI models the ionic states of 
each abundant element (e.g., C, N, O, Ne, Mg, Si, S, Ca, 
Fe, Ni) are solved for through time-dep endent ioniza- 
tion equations. APEC (Smi th et al.ll200lf ). which is used 
throughout this work, on the other hand does an Ioniza- 
tion Equilibrium calculation for X-ray emission spectra 
of collisionally-ionized diffu se gas calculated usin g the 
ATOMDB code v2.0.1 ( |http: / /www.atomdb.org/[ ) that 
includes a set of relevant ions and their lines and per- 
form a population calculations. The atomic data and 
line transitions information are at a more comprehensive 
state than other models like Mekal and has a wavelength 
coverage from IR to X-ray. 

However, when the resulting grid of feasible values 
for To is sampled it was noted that the data does not 
constrain it at all. In case of SN 2004dj, the post 
reverse shock density at the first observation epoch is 
~5x 10 6 cm -3 , and the age is 5 x 10 6 s. This will give 
a To > 10 13 s cm~ 3 , which should be enough time for 
the plasma to come to ionization equilibrium. This may 
not be strictly true for the material behind the forward 
shock, however it would be too hot and too low den- 
sity to contribute significantly to the Chandra spectrum. 
Therefore it is safe to make the assumption of ionization 
equilibrium while modeling the Chandra X-ray spectrum 
of SN 2004dj , although higher resolution data such as will 
be available with Astro-H might reveal NEI features. 

5.3. Particles and Magnetic Fields 

Radio or X-ray light curves of supcrnovae are often ex- 
plained using varying amounts of accelerated electrons 
and amplified magnetic fields. However it is difficult to 
ascertain the relative contribution of each. By modeling 
the observed X-ray spectrum as a combination of thermal 
and inverse Compton flux, we have been able to deter- 
mine the equipartition factor for SN 2004dj. Our result 
shows that the plasma is away from equipartition, with 
electrons having more of the total energy than magnetic 
fields. Since e e and cb have also been found individually, 
theories of particle acceleration and magnetic field am- 
plification can now be tested against them. This is an 
unique opportunity. 

Since the electrons and magnetic fields are responsible 
for producing the synchrotron spectra of supernovae, it 
is important to lo ok breaks in th e radio spectrum of a 
young supernova. iChandra et al.l (j2004f ) associate such 
a break in SN 1993J (not a type IIP, but a type IIB su- 
pernova) with the phenomenon of synchrotron aging of 
radiating electrons. The magnetic field in the shocked 
region can be calculated from the break in the spectrum, 
independent of the equipartition assumption between en- 
ergy density of relativistic particles and magnetic energy 
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density. The synchrotron cooling break appears at a cer- 
tain frequency when electrons radiating above that fre- 
quency have lost a significant portion of their energy to 
synchrotron emission. Demanding that the synchrotron 
loss time scale is equa l to the age of the sup ernova, we 
can use Equation 5 of iChevalier et ail (|2006| ) to express 
the frequency of the synchrotron cooling break as 
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which has the sam e temp oral evolution as derived in 
iChakraborti fc Ravi (|2011l) . For the epochs at which 
it was observed in the radio This is much higher than 
the frequencies accessible with the EVLA, therefore syn- 
chrotron cooling will be unimportant for studying the ra- 
dio light curves of SN 2004dj. However for future studies 
of young type IIP supernovae, these frequencies may be- 
come detectable with mm wave observatories such as the 
SMA or ALMA. 

Another important mechanism for electron cooling is 
inverse Compton losses against the low energy photons 
from the supernova photosphere. Due to this process 
high energy electrons may lose bulk of their energy in 
boosting the seed photons to X-rays, thereby contribut- 
ing to the X-ray flux, as discussed in this paper for SN 
2004dj. Demanding that the inverse Compton loss time 
scale is equal to the age of the sup ernova, we can use 
Equation 6 of IChevalier et ahl (|2006T ) to express the fre- 
quency of the inverse Compton cooling break as 
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10 4 kms-V V 10 42 ergs s~ 1 J K 1 

For the determined relevant values of the parameters de- 
termined in this work, the v\c is within the frequency 
coverage of the VLA. Therefore, any interpretation of 
the radio light curve of SN 2004dj should account for in- 
verse Compton losses and expect to see such a cooling 



6. CONCLUSIONS 



In this work we have used SN 2004dj as a prototype for 
type IIP supernovae to answer some fundamental ques- 
tions about their circumstellar interaction. 

1. What is the origin of X-rays detected by the Chan- 
dra from young type IIP supernovae? The X-rays 
detected by Chandra arise- from a combination of 
thermal and non-thermal processes. They are dom- 
inated at early times by optical supernova photons 
inverse Compton scattered by relativistic electrons 
at the forward shock which falls off as oc L^/t. 
At late times it is dominated by thermal emission 
from reverse shock heated plasma which only falls 
as oc t~ l , 

2. How fast was the progenitor losing mass via winds 
in the final phase before explosion? The progenitor 
was losing mass at a rate of M = (0.32 ± 0.11) x 
10~ 6 (Vw/10 km s" 1 ) M Q yr" 1 . This is less than 
the value expected from mass loss prescriptions in 
the literature for the putative 15M0 progenitor. 

3. How efficiently does the supernova accelerate cos- 
mic ray electrons in its forward shock? Around a 
third (e e ~ 0.39) of the energy thermalized by the 
collision of the ejecta with the circumstellar mat- 
ter is used in accelerating electrons to relativistic 
velocities. 

4. What is the extent of turbulent magnetic field am- 
plification in the post-shock material? Around a 
tenth (ee ~ 0.082) of the thermal energy available 
is used in the turbulent amplification of magnetic 
fields. 

5. Are synchrotron and inverse Compton losses im- 
portant in explaining the radio lightcurves? In 
the case of SN 2004dj synchrotron losses are unim- 
portant while inverse Compton losses will have to 
be taken into account while explaining its radio 
lightcurve. 



This research has made use of data obtained from the 
Chandra Data Archive and software provided by the 
Chandra X-ray Center (CXC) in the application pack- 
ages CIAO and ChlPS. 
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